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Quark-Model Baryon-Baryon Interaction Applied to the 
Neutron-Deuteron Scattering (II) 

Polalization Observables of the Elastic Scattering 

Kenji Fukukawa and Yoshikazu FuJiWARA 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

The neutron-deuteron (nd) scattering is solved in the Faddeev formalism, employing the 
energy-independent version of the quark-model baryon-baryon interaction fss2. The differen- 
tial cross sections and the spin polarization of the elastic scattering up to the neutron incident 
energy En ~ 65 MeV are well reproduced without reinforcing fss2 with the three-body force. 
' The vector analyzing-power of the neutron, Ay (6), in the energy region E„ < 25 MeV is 

largely improved in comparison with the predictions by the meson-exchange potentials, thus 
, yielding a partial solution of the long-standing Ay-puzzle owing to the nonlocality of the 

. short-range repulsion produced by the quark-model baryon-baryon interaction. The large 

' Coulomb effect in the vector and tensor analyzing-powers in E„ < 10 MeV is also analyzed 

based on the Vincent and Phatak method and recent detailed studies by other authors. 

in : 
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§1. Introduction 

O ' The QCD-inspired spin-flavor SUq quark model (QM) for the baryon-baryon 
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interaction, devefoped by the Kyoto-Niigata group, has achieved accurate descrip- 
tion of available nucleon-nucleon (NN) data, comparable with the modern meson- 
exchange potentials.!^ The naive three-quark structure of the nucleon is incorporated 
^ ■ in the microscopic framework of the resonating-group method (RGM) for two three- 

i quark clusters, leading to the well-defined nonlocality and the energy dependence of 

the A^A^ interaction, inherent to the RGM framework. In particular, the short-range 
^Nji' . repulsion of the A'^A'^ interaction is mainly described by the quark-exchange kernel 

originating from the color-magnetic term of the quark-quark interaction. The energy 
dependence of the interaction is eliminated by the standard off-shell transformation, 
utilizing a square root of the normalization kernel.--' This procedure yields an ex- 
tra source of nonlocality, whose effect was examined in detail for the three-nucleon 
(3A^) bound state and for the hypertriton.® The QM baryon-baryon interaction thus 
^ . constructed is expected to give quite different off-shell properties from the standard 

^ I meson-exchange potentials. It is therefore interesting to examine predictions by the 

QM A'^A'^ interaction to the 3A^ scattering, especially in this renormalized framework 
with no explicit energy dependence. 

In a previous paper,® referred to as I hereafter, we have developed a new al- 
gorithm to solve the Alt-Grassberger-Sandhas (AGS) equation^ with the deuteron 
singularity, employing the Noyes-Kowalski method.'^''^ Another notorious moving 
singularity of the free three-body Green function is treated by the standard spline 
interpolation technique developed by the Bochum-Krakow group.'^"'^ The AGS 
equation is solved in the momentum representation, using the off-shell RGM t- matrix 
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obtained from the energy-independent renormalized RGM kernel. The total cross 
sections derived from the optical theorem and nd elastic differential cross sections 
are well reproduced up to the neutron incident energy En = 65 MeV in the labora- 
tory system. It was found that the predicted elastic differential cross sections have 
larger diffraction minima than the experiment in the energy region En = 35-65 
MeV, in contrast with the predictions by the meson-exchange potentials. This is 
in consistent with the fact that fss2 reproduces nearly correct triton binding en- 
ergy without the three-body force. In these calculations, we have used the Gaussian 
nonlocal potential constructed from the original fss2 interaction.'^^ This potential 
preserves not only the on-shell properties of the t-matrix with the accuracy of less 
than 0.1 degree for the phase shift parameters, but also the nonlocality of the inter- 
action completely, which was confirmed by recalculating the triton binding energy 
with this potential.^' The calculation of the S'-wave nd scattering length by this 
potential yields almost correct values ^Ond = 0.66 fm for the spin doublet scattering 
length and ^a„d = 6.30 fm for the quartet scattering length, without the three-body 
force and the charge dependence of the A^A^ force.^^ 

In this paper, we extend the study of the nd elastic scattering in I and Ref. llSp 
to various types of spin polarization observables, and compare them with the exper- 
imental data and predictions by other theoretical calculations. In the low-energy nd 
and pd scattering, there is a long-standing nucleon analyzing-power puzzle, which 
implies the failure of rigorous 3A^ calculations to account for the magnitude of the 
measured analyzing-power 

j^^^ff^mm We find that there is no serious discrepancy 
in Ay (6) at the low energies En < 25 MeV, although the difference of about 15% still 
remains for the magnitude of maximum peaks. We have analyzed the low-energy nd 
and pd eigenphase shifts below the deuteron breakup threshold and found that the 
main origin of this improvement is the more attractive feature of the phase 
shift in our model, originating from the nonlocal description of the short-range repul- 
sion.fl^ The comparison of tensor-type deuteron analyzing-powers with experiment 
is blurred with the Coulomb effect, since all the experimental data are for the pd 
or dp scattering. We here apply the Vincent and Phatak methocP®''^ using the 
cut-off Coulomb potential to the pd scattering, and show some preliminary results 
in comparison with recent detailed studies by other authors in the Ep < 10 MeV 
region.'^'t^ We find that this method works well to reproduce the characteristic 
behavior of the forward angular distributions for the dp tensor analyzing-powers. 

The organization of this paper is as follows. In §2.1, we first develop a general 
framework to calculate the spin polarization observables for the nd elastic scattering, 
using the Blatt and Biedenharn technique.'23 The applications to the analyzing- 
power and polarization of the nucleon and deuteron, the polarization transfer for 
the nucleon, the nucleon to deuteron polarization transfer, and the spin correlation 
coefficients are made in subsections 2.2, 2.3, 2.4, and 2.5, respectively. The definition 
of the deuteron spin operators, used in this paper, are summarized in Appendix A, 
together with the relationship between various different notations. The results of 
polarization observables are shown in §3. In §3.1, we discuss the nucleon analyzing- 
power, focusing on the analyzing-power puzzle. The vector and tensor analyzing- 
powers of the deuteron are discussed in §3.2, where the Coulomb effect is important at 
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Ep < 10 MeV. In § 3.3, we extend the analysis of the low-energy observables below 
the deuteron breakup threshold, developed in Ref. [T3|) . by employing the present 
approach to the Coulomb force. Various types of polarization transfer and spin 
correlation coefficients are discussed in § 3.4 and § 3.5, respectively. The last section 
is devoted to a summary. 

§2. Formulation 

2.1. General framework 

In this subsection, we derive general formulas for the spin polarization observ- 
ables for the nd elastic scattering by using the Blatt and Biedenharn technique.'^ 
This method is used in the previous paper I to calculate the differential cross sec- 
tions and the extension to the various spin observables is rather straightforward. We 
calculate I = Tr {fSif^Sj}, where / is the nd scattering amplitude defined by the 
partial- wave amplitudes f^^is') [iSc) Eq- 1(2.92^3 as 

/(Q/, Qi) = 47r Yl fi'S'Mes.) Yj ^i^'S'^VJ. (9/; spin) yits,)jj^ (9*; spin) , 

(2-1) 

where 3^(f5^)jj^(q; spin) = [Y£(q)x5c (12; 3)] jj^ is the angular-spin wave function. 
The spin wave function is defined by XScSczCl^; 3) = [xi(12)xi (3)]s'cSc3- The spin 
operators Si and Sf are in general expressed as a tensor operator defined by 

Jz 

where s = or 1 and A = 0, 1, or 2. The neutron spin operators {l,cr} and the 
deuteron spin operators {1, S', S"^^^} form 4 x 9 = 36 independent matrix components 
in the spin space spanned by XScScz (12; 3) with 5c = 1/2 and 3/2. The deuteron spin 
operators S and S^'^'^ are summarized in Appendix A, together with the relationship 
between various notations. The standard expression for the nd scattering is in the 
Cartesian representation such as aaS^^, which are expressed by the linear combina- 
tion of the tensor coupled form in Eq. ()2-2p . To calculate / = Tr {fStpSf}, we use 
Eq. ()2-ip and take the trace for the spin variables. Then, we obtain 

I = Tr{fS,pSf} 

= (4vr)2 Y Y ^i'S'J{iSc){f(e^s^){ISc)} 

{i'S',)iis,)j(^e'S',)iISc)J 

X Y 0^{i^Sc)JJz (Q*; spin)|5i|3^(^^) jj^ (q^; spin))spin 

JzJz 

^^F'5^)JX(^/;spin)|5/|3^(£/s/)jj^(gj;spin))spin , (2-3) 

*-* In the following, we cite equations of the previous paper I with adding I in front of the 
equation number. 



cr 



(s) g{X) 



(2-2) 
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where the last two matrix elements are taken only for the spin variables. We assume 
that the spin operators Si and Sf are tensor operators in Eq. ()2-2p and calculate 



/(Vi;A//i/)=Tr{/5(^)/t5g^^} . 

We first calculate 

Jxf^iQt) = (3^(^Sc)JJ.(9i;spin)|5i-^^|3^(j§^)jj^(q,;spin))spin , 
in the standard recoupling technique. It is rather easy to derive 



(2-4) 
(2-5) 



LM 

xZj^\£jIj;SAL) , 



(2-6) 



where the recoupling coefficient {iJiJ; ScScL) is defined by 
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Z\^\UU; ScScL) = {-Y+J-S- i J {eoiO\LO) 
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£ Sc 


J 





(2-7) 

using the unitary form of the 9-j coefficient and the unconventional reduced matrix 
element (xsdl'^^'^'^llx^ )unc- We can use this formula both for 5^^'^ and S^/^ in Eq. 



()2-3p . changing toqj etc. Then the sum over and Jz can be taken. We further 
use the symmetry property of z\^\lJlJ; ScScL) in Eq. (j2-7p : 



Zl^\ej£J;ScScL) = {-)^+''+^+^^-^^ Zi^\£JU; ScScL) , 
by assuming 5^'**'*^ = {—)'^S^^. After all, we obtain for Eq. I^i 



I{Xifii;Xffif) = ( 

ie'S',)i£S,)j(^i>s',){ISc)J 



(2- 



fie'S',)ieSc){f(e's',){Is,)} 



LM,L'M' Kv 

X Z^^^\uIJ; ScScL) Zi^f^ {I'Ji'J; S'^S'^L') . (2-9) 



We write this expression as 

I{Xi^l^■,Xf^If) = Y Yl (")^'~^' A'^'soc^s.) {/(i)50(F5e)} 

X Y Ci:^''\h^^^Af^^f.qf,qi) z\t^\UU-ScScL) Z^^^\i'jI'J-S^ (2-10) 



L,L',K 
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The spatial part C^'^ {Xi^i, XfiJf;qj,qj} can be written in various forms: 

MM'u ^ ^ 



(-) 



Xi+Xf+L+D 



but the separation into the invariant part and the non-invariant part is not easy 
except for the / = and 1 cases. We should note that the L and L' sum in Eq. (|2-1U|) 
is only for \L — L'\ =even because of the parity conservation. We use a simplified 
notation C^' {Xifii, Xj^f) = C^' {Xi^i, Xffif;qj,qi) in the following and choose a 
special coordinate system with = (^,0) and = = (0,0), where the z-axis is 
the beam direction and 6 is the scattering angle in the center-of-mass (cm) system. 
We use 

Yemi^z) = ^^m(0,0) = 6m,0 



47r 



m' + lm'l f (i'-\m'\)l l^'l 

y„„,(«,0) = (-)^^^L_L_^ «) , (242) 

which yield a basic symmetry property 

Y{e'e)KuiQf,ez) = (-)"+'' Y^(,,iy_^{qf,e,) 

= (_)^'^(^W-z.KO)y|^^P,';(cos0) for ^>0. (2-13) 
If we use Eq. (|2-12p in the first expression of Eq. ()2-lip . we find 

c^^'(A,^„A/^/) = {-)^^+>^^+^f+>^fct'''{Xi,-^,i,Xf,-^,f) 



= i-rf{Kti,x, - /i,|Lo)(K/x,A/^/|LV. + p{;;+^^(cos e) 

for Hi + fif>0. (2-14) 
If we set here Xffif = 00, then we obtain 

C7i"^'(Ai/.,,00) = (-)^'+A'»C^^'(Ai,-/.,,00) 



Sn,L'{L'^I^Xi - ^,|LO)^/|^J-p^ P^;(cos e) for /i, > . (2-15) 



For Xifii = 00, we find 

ct'^'mXffif) = (-)A/+/^/c^^^>M/,-^/) 

= 6^,Lj{L'fifXf - l^fimJ fj^, ~ ^^j; PlHcos 9) for > . (2-16) 
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If we further set Aj/ij = 00, we gain 

^^■^'(00,00) = S^^lSl,l' Pl{cos 6) , (2-17) 
which yields the differential cross sections in Eq. 1(3.5) by using the reduction 

) {UU; ScScL) = 5k,l5s^ Z{£jIJ; ScL) . (2-18) 
The reduced matrix element of S^^^ in Eq. ()2-7p is given by 



ixsj 



is) g{X) 



a 



^s+X-f 



if) 



A s f 



1 ^ Sc 



with (|||cr||i),nc = V3, (1||S||1) 



V2, and (l||5(2)||l)un 



2 / unc ) 



(2-19) 
5/3 (see Eq. 



2.2. Analyzing-power and polarization 

The analyzing-power is characterized by A//i/ = 00 and the spatial function in 
Eq. ()2T5p is used. The nucleon vector analyzing-power is defined with Si = cr. The 
expression in Eq. (|2-15p with Aj = 1 



C,^'^'(lM,00) = (-)i+A^C,^.^'(l,-;.,00) 
= 6^^L'{L'fil- ^l\LO)J^^fy^^ P^,{cos9) for fi > 



(2-20) 



is transformed to the Cartesian representation by the spherical vector rule in Eq. 
(|A-3p . If we write these as C^' (a), we find that x and z components are zero 
because of the parity symmetry in Eq. (j2-20p and (L010|L0) = 0, respectively. (Note 
that we only have L' = L for Aj = 1.) The y-component is from (Lll — l|i^O) = l/\/2 
with L > 1 



Pi, {cos 



(2-21) 



VL'{L' + 1) 

The nucleon analyzing-power Ay{9) is usually defined through 

where n = [q^ x Q/]- Since n = ey, we find that Ay (9) is calculated from ly/Io with 
Eq. ()2-2ip . where Iq = Tr {//^}. The reduced matrix element for Z^j^\iJ£J; ScScL) 



(2-22) 



IS 



s. 



- I 1 

I Sc 1 



(2-23) 
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The vector analyzing-power of the deuteron is similarly calculated, by setting 
Si = S. Only the reduced matrix element is different, which is explicitly given by 

{xsJ\S\\XsJunc = i-lf^^^V6^c!^^l ~ }. (2-24) 

From the definition of iTu in Eq. ()A-22p . we should multiply the factor \/3/2, namely, 
iTu = (\/3/2)(/(y)/Io). 

The tensor analyzing-power of the deuteron is obtained from the Aj/Zj = 2fi 
case in Eq. ()2-15p . From the parity conservation embodied in Eq. ()2-15p . the two 
components with 2// and 2, —fj, give the same contributions. Thus the factor \/3 
in Eq. (|A-16p yields = \/3(/(2/i)//o), where the reduced matrix element for 
Z^}{£J£J; ScScL) (in this case, L' = L, L zt 2) is given by 

(xsJ|5(2)||xg>nc = (-l)^^+^V5le{ J (2-25) 

The spatial functions /(2/i) = /(2/i, 00) are calculated from 
C^^'(2/i,00) = i-r C t'''{2,-fs, 00) 

= 6^,L' {L'fi 2, -/i|LO)^|^J^^ P^,(cos 9) for ^ = 0, 1, 2 . (2-26) 

It is convenient to write the final results for the analyzing-power in the form 
similar to the differential cross sections in Eq. 1(3.5) (We take the sum over 5"^ and 
Sc inside.) 

Io = Tr{fp} = 6ao = Y,BL Piicose) , 

L 



I{2m)=Tr{fS(^^p} = J2B^^^„ 



\l' - m)\ 



(L' + m)! 



x(L'm 2,-m|L0) P^(cos 6*) for m = 0, 1, 2 , (2-27) 

where P™(x) (m = 0, 1, 2, • • • ) are the Legendre by-polynomials and the II sum in 
I {2m) is only for L' = L, L ±2 from the parity conservation. The coefficients, B^, 

B^^^ and B^^^^,, are given by 
Bl' ^ 

^(2)^ I ii's^)ies,)j(^i,s',){ISc)J 

L.L ' 



X < 
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6g^^^Z{eJiJ;ScL') 

i Z^^\mj; ScScL') } Sg, g,Z{£'J£'J;S',L') 
zf;}{UU;ScScL) 



(2-28) 



The explicit expressions of the CG coefficients aUow us to express /(2m) (m 
0, 1, 2) in Eq. (f2^ as 



/(2m) = Tr 



/•} = E 



L>m 



1 



Pr(cos e) , (2-29) 



I y/(L-l)L(L+l)(L+2) 



for m = 0, 1, 2 with 



L(L + 1) 



(2L-l)(2L + 3) 



4^2 + 



/ 3(L + l)(L + 2) (2) 
2(2L + l)(2L + 3) ^+2'^ 



5i 



+1 



3(L - 1)L ^(2) 



2(2L-1)(2L + 1) ^-2,L ' 



' 3 


2(2L- 


l)(2L + 3) 


liL- 


-1)(L + 1) 


V(2L- 


-1)(2L + 1) 


' 3(L- 


l)(L + 2) 


2(2L - 


l)(2L + 3) 



4^2 + 



L{L + 2) 



(2) 



(2L + l)(2L + 3) ^+2'^ 



5 



(2) 

L-2,L ' 



(L-l)L 



+ 2 V (2L + 1)(2L + 3) ^^+2'^ 



5 



(2) 



1 / (L + l)(L + 2) (2) 
2V (2L-1)(2L + 1) ^-2.^ 



(2-30) 



The polarization of the outgoing particles is calculated from the Aj/ij = 00 case in 
Eq. ()2-16p . In the cm system, the vector polarization of the nucleon or the deuteron 
is the same as the analyzing-power. Namely, we find (00, y) = C^'"^ (y,00). For 
the tensor polarization of the deuteron, we find 



C^^'(00,2/x) = i^C^^'(2^,00). 

1j 



(2-31) 



2.3. Polarization transfer for the nucleon 

This is the case with Si = a and Sf = cr. Setting Aj = 1 and A/ = 1 in Eq. 
(l2Ta . we find 
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for Hi + Hf>0. (2-32) 

From the symmetry of the parity conservation, the independent types of the polariza- 
tion transfer are very restricted. We use a simphfied notation C^,;^' = C^' (1^, 
and express this symmetry as C^,^' = (—1)^+'^ C-^^-^i from Eq. ()2-32p . For the fixed 
values ^ii + fif = fz = 0, 1, 2, we have the foUowing five combinations of fii and 
fif. = (0, 0), (1,-1) for /, = 0, (10), (01) for = 1, and (1,1) for = 2. 

The correspondence to the Cartesian representation is found, for example, as 

Kx,x ~ CTx^x = 2^~^l + 0'-l)(-<7l + CT-l) ~ 2^'-^^^ ~ ^I'^i ~ C*-!,! + C-i,-i) 
= 2^^^^ ~ ^I'-i ~ Ci^^i + Ci^i) = Cii — Ci,_i , 

= , 

~ Cra:<7^ = — ^(-(Ti + Cr_i)(To ~ -^(-Cio + C-i^o) = "V^ Cio , 

(2-33) 

We obtain 

Kz,z ~ Coo ; -f^z.x ~ Cii — Ci_-_i , Ky^y ~ — (Cii + Ci _i) , 
Kx,z ~ (-V2) Cio , Kz,x ~ (-\^) Coi , 

Kx,y = Ky^x = Ky^z = Kz,y = . (2-34) 

After all, five independent polarization transfers of the nucleon are calculated from 
the following spatial integrals: 

Kz,z ~ Ci"-^'(10, 10) = (k010|L0)(k010|L'0) Pl'(cos 9) , 

Kx,z ~ (-^/2)Ci"^'(ll, 10) = -V2{k 1 1, -l|LO)(KllO|L'l)^=i==P2,(cos 6) , 

^/ L [L +1) 

Kz,x ~ (-^/2)C,^'^'(10,11) = ^/2(At010|L0)(A^011|L'l)^==i== p2,(cos 9) , 

^/ L [L +1) 

Kx,x \ f C^'^'(11,11) -Ci^'^'(ll;l,-1) 
-Ci"^'(ll,ll)-Ci"^'(ll;l,-1) 



K 



= (k11,-1|L0) ^{t^in\L'2)^^^y^Pl,{cos 0) + (k 1 1, -1|L'0)Pl'(cos 6) 
The reduced matrix elements of the two Z-factors are for a, given in Eq. (|2-23|) 



(2-35) 
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We should note that the polarization transfers in Eq. ()2-35p are for the cm 
system. These are however usually defined in the lab system and the notation K^^ 
with the upper /3' is used to specify the /3-axis in the laboratory system. We have 
to rotate -f^a,^ with respect to the coordinate (3 for the outgoing particle. In the 
laboratory system, the scattering angle of the nucleon, 9 in Eq. (j2-12p is transformed 
into 9i, which can be calculated from 

sin 9 

t8Ln9i = — —. (2-36) 

W 



7 cos u + 



Here, 7 and (3 are the usual relativistic factors for the cm to lab transformation, and 
/3* is the /3 factor of the outgoing nucleon in the cm system. For the neutron-incident 
nd elastic scattering these are calculated from the "kinetic" energy of the incident 
nucleon Tiab = Eiab — M^c^. {M]\f is the nucleon mass.) If we assume the masses of 
the nucleon and deuteron as mi = and m2 = 2Mj\f, the necessary factors in Eq. 
(j2-36p are calculated to be 



with X = Tiah/{M]\jc^). The rotation of the Cartesian frame (x, y, z) in the Madison 
convention around the y-axis by Oi yields the transformation of the unit vectors 

e'^ = Sz cos 9i + Ex sin 9i , 
e'j, = —Sz sin 9i + cos 9i , 

e'y = By . (2-38) 

The components of the spin vector a of the outgoing nucleon are therefore trans- 
formed to 

(^z' = COS 9i + ax sin 9i , 
Cx' = — o"z sin 9 1 + ax cos 9i , 

ay' = ay , (2-39) 

since a a' = {cr ■ e'„) and iTq = (<t • Cq). In order to write down the explicit expression 
for Ka , we note that K^^fs in the cm system is actually functions of pj = {9,(p) 
under the assumption of = = (0,0). We therefore understand Ka,/3 in Eq. 
(|2-35p are actually K^^p = Ka^(s{6,0), using the notation Ka^i3{9,ip). Taking these 
into consideration, we obtain the following results. 



Kx' 
Kx' 



Kz,z COS 6*1 + Kz^x sin 9i , 
-Kz^z sin 9i + Kz,x cos 9i , 
Kx,z COS 6*1 + Kx,x sin 9i , 
-Kx. sin 9i + Kx,x cos 9i , 



Kyy = Ky,y , (240) 
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where -K'a,/? = Ka^j3{9, 0) and 9i is calculated from Eqs. ()2-36p and ()2-37p . 

In fact, we further need smah relativistic corrections for the spin directions, 
as discussed in the A^A^ spin transfer coefficients. (See for example Ref. l28p and 
references therein.) We however neglect these, since we deal with the low-energy 
scattering for the time being. 

2.4. Nucleon to deuteron polarization transfer 

The vector-type nucleon to deuteron polarization transfers in the cm system 
are calculated from the same spatial functions as in Eq. (j2-35p . but with a different 
reduced matrix element of S given in Eq. ()2-24p for the second rank-one Z factor. The 
same notation is used for these polarization transfer coefficients with the vector- 
type deuteron polarization. However, we should be careful with some complications 
originating from the fact that we are now detecting the recoil particle. The scattering 
angle 62 of the recoil deuteron is given by 



namely, O2 = tt /2 - Q with 7r/2 corresponding to ^ = and to = vr. The direction 
of 62 is opposite to ^i, so that the y' axis is not equal to the original y axis but to 
the opposite to it. (e^ = — e^). If the nucleon comes out on the left-hand side of 
the beam direction, the deuteron turns to the right-hand side. This changes the sign 
of e'y. This definition seems to be inconvenient to the experimentalists, since the 
outgoing deuteron is further turned around by the spectrometer magnet to measure 
the polarization. To circumvent this difficulty, we rotate the whole system around the 
original z-axis by vr. Then the y-direction becomes the original By and the deuteron 
comes out to the left-hand side (the same direction as the nucleon before). We can 
use 



for the deuteron spin operator S, similarly to Eq. (|2-39p . but we should be careful 
to use Ka^p{9,T:) and not Ka^j3{9,Qi). We therefore find 




(2-41) 




Sz cos 92 + Sx sin 92 , 
—Sz sin 6*2 + Sx cos 6*2 



(2-42) 



Kz' = Kz,z(e,7:) cos 92 + i^.,x(0,vr) sin 02 , 
i^/' = -Kz,z{G.^) sin 92 + Kz,x{e,7:) cos 92 , 
Kx'^' = Kx,z{0,tt) cos 6*2 + Kx,xiO,TT) sin 6*2 , 
Kx""' = -Kx,z{0,Tr) sin 6^2 + Kx,x{0,ir) cos 6*2 , 
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See Sec. 5-3-2 in page 751. In our expression, this symmetry appears in Eq. ()2-34p 
as the fact that the coefficient C^,,^^ changes the sign when /ij + //^ = odd. This is 
because the Yii^'iO,'^) in Eq. (|2-12p gives an extra factor (— l)^'i+^'/ from the wave 
function <Pjn'{'p) = e"^"^ / V2tt . For the polarization transfer coefficients this phase 
factor is the same as the even-odd character of + Ny, where A^^^ and Ny are the 
number of times x and y appear. Namely, from Eq. ()2-34p . we find that only K^^z and 
Kz^x change the sign, since fj,i + fif = odd. In other words, Kx^z{(), tt) = —Kx^z{S, 0) 
and Kz^xid,'^) = —Kz,xiS,0), the others no phase change. If we use the simplified 
notation i^o,/3 = Ka^p{6,Q) as in Eq. (|2-40p . we finally obtain 



l\y 



Kz,z COS 6*2 - Kz,x sin 6*2 , 
-Kz^z sin 6*2 - Kz^x cos 6*2 , 
-Kx,z COS 02 + Kx^x sin 6*2 , 
: Kx^z sin 02 + Kx,x cos 02 , 

Ky,y . (2-44) 



Note the phase change in addition to 0i -^02- 

For the tensor-type deuteron polarization, we consider the matrix elements for 

(2) 

Si = (7^. and Sf = S^i^ . Since Aj = 1 and = 2, we deal with the spatial integrals 
of the type 

Ct'^'il^ii, 2^;) = (-)i+^»+^/C,^'^'(l,-M,; 2,-^;) 



= (-r/(K/xa-/.i|L0)(KM.2^;|LVi + /./)^ |^i^|]i^|^^ji Pi:;+^^(cos 0) 

for ^i + ^j = 0, 1, 2, 3 . (2-45) 

We again use the notation C^,|j' = C^'^ (1^, 2fi') and express the parity symmetry as 
Cf_i,f_i' = (— 1)^^'^"'"^'C_^___^/ fromEq. ()2-45p . From the parity conservation = 
1, we only have L' = L, Lib2. We can calculate independent non-zero coefficients in a 
way similar to the nucleon polarization transfer. Namely, we express the tensor-type 
nucleon to deuteron polarization transfer coefficients Ka"^ from the spatial matrix 
elements for the corresponding spin operator aaVp^^y. The transformation from 7^/3,7 

('2') 

to si:> is given by Eqs. ()AT9p and ()AT3p . Again, from parity conservation, some of 
the coefficients are zero. In particular, we have a relationship Ka^ + Ka^ + Ka^ = 
for all a = X, y, z, since Vx,x + 'Py,y + 'Pz,z = 0. Because of this symmetry, only 
seven coefficients are independent. These are 

i^,-'^~Z-|(Ci2-C_i,2), 



-V3C,o + ^C,2 + ^C_,2) 
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Ry'y ~ i 



V3Cio - ^Cu - 



K^^' = - {K^^^ + K^y) ~ i 2^/3 Cio , 



K^'y^ii-3)Co2 , 



(2-46) 



The correspondence also yields 



Ry^y 



00 



Rl 



-2)R: 



x,x 



V6Q 



00 



(2-47) 



However, these are all zero since Cqo = {—Y^^'^^Coo = from the parity conserva- 
tion. The independent seven spatial functions are the linear combinations of Ci-i 
for = 0, Cio, Coi, C_i,2 for = 1, C02, Cn for = 2, and Cu for = 3. Af- 
ter all, we calculate seven independent tensor-type nucleon to deuteron polarization 
transfer coefficients by the spatial integrals 



Ky'"" -i — 

V2 



(11,22)-C7^^ (1,-1,22) 



C,^'^'(ll,21)+C^^'(l,l, 2,-1: 



Ry' ~i 



Ry^y ~ i 



y 



- V3C,^'^'(11, 20) + ^^,^'^'(11, 22) + -1, 22) 

-^/3C,^'^'(11,20) - -|c,^,^'(ii,22) - -|cf,^'(i,_i,22) 
(i^^-^ + JC^'J') ~i2V3C^'^'(ll,20) , 



^ V2 L 



Cf"^'(ll,21)-C^'^'(l,l, 2,-1) 



R'fy ~ i(-3)C7^'^'(10,22) 
~z3C^'^'(10,21) . 



(2-48) 



The reduced matrix element for the second rank- two Z- factor is given by Eq. ()2-25p . 

The transformation to the laboratory system can be carried out in a way similar 
to RJ'. We first use Vaa = 3S'q,^ — 2 and Eq. (|A-15p . and derive spin rotation rule 
according to Eq. ()2-42p (we here use the simplified notation 9 for 62 in Eq. (|2-42p .): 

= ^x.x.(cos Of + V,,{sm ef - P^,(sin 20) , 

V' — V 

' yy — ' yyi 

T^L = 7^xx(sin Of + P^^(cos Of + P^,(sin 26) , 



'^xy ~ "^xy COS 



Vyz sin 
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I yz 



Vxy sin 



Vyz cos 



Kz = i'Pxx - Vzz) 2 (sin 26) + V^cos 29) 



(2-49) 



The vr-rotation around the z-axis yields 



K ^ 


,y 


= K^x'^if' 


,^)cos 02 -i^^'^(^,^)sin 02 , 




k/ 


,z' 


= K'HO 


,^)sin 02 + Ky''{e,TT)cos 02 , 




k/ 


,y' 


= K'fy{6 


,^)cos 02- Ky''{e,TT)sm 02 , 




k/ 


,z' 


= K^'y{e 


,^)sin 02 + Ky''{e,TT)cos 02 , 






,z' 




9,7t)-K^^'{0,7t)] ^(sin 2^2) +i^y"'^(0,vr)(cos 2^2) 




I\y 


,x' 


= K^^^{e 


,^)(cos 02)^ + K^''{e,TT){sm 02)^ - K^''i0,TT){sm 


2^2 


Kyy' 


,y' 


= Ky'y{e 


,71") , 




Ky'' 


,z' 


= K^^'-'ie 


, ^)(sin 62)^ + K^^'{0, 7r)(cos ^2)^ + K^''{0, vr)(sin 


2^2 



(2-50) 

If we further use the odd-even character of Ka^{9, 0), we find that Kx'^ , Ky'^ , -^y'^, 
Ky'^ , and Kz'^ (namely, A^^ + Ny = odd) change the sign from Eq. ()2-46p . We 
therefore obtain the following final result. 



K y 



Kz 
Kz"^^ 



K ^ 

J\y 

K 

l\y 

K y 

JXy 

K.„ 



■ -K^'y COS 02 - Ry^' sin 02 , 
-Kl^y sin 02 + Kl'' cos 02 , 
i^f ^ cos 02 + Kl'' sin 02 , 
i^f^sin 02-Kf'cos 02 , 

[-K-y^- + Kl^'] i(sin 202) + K,^'^(cos 2^2) , 

: -i^;^'"(cos 02? - K^'^isin 02)^ - K^^'isin 202) , 

_ Ky'y 

^^y ' 

-K;^'"(sin 02)' - K^"'"(cos 02)' + K;^'"(sin 2^2) , 



(2-51) 



where K^'^ = K^''^{0,O) and 02 is calculated from Eqs. (|2^ and (12^71) . Note that 
Ky'^ changes the sign. 



2.5. Spin correlation coefficients 

In this case. Si = [o'S]^^^ or 5, = [crS^^^jli^^ together with 5/ = 1. The exchanged 
case, Si o 5/ is considered in a similar way (except for the phase correction of the 
spin reduced matrix elements related to the symmetry in Eq. (|2-8|) ). We define the 
matrix element /(A^), which is a special case of Eq. ()2T0p with Aj/ij = A/i and 
Xffif = 00: 

/(A/z) = I(A/x;00) = rr{/5W/t} = (-)A+mj(a, 
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= XI Yl ^'^^"^ fi'S',){iSc){f^I'S',)(Is,)} 

X ^C7^;^'(A^,00) z[^\£jIJ; SjcL) 6^, Z{e' Ji'J; S'^L') . (2-52) 

L,U 

Here, the spatial functions C^l^ (A/u, 00) are given in Eq. ()2-17p for A// = 00, in Eq. 
(|2-20|) for A = 1, and in Eq. (|2-26|) for A = 2. The reduced matrix element for the 
first Z^j^^ function is given by Eq. ()2-19p . 

Let us first consider 5^^^ = [ct^jIj^^ case. The reduced matrix element from Eq. 
(|2T9]) is 



(X5j|[cx5]W||X5>nc = (-l)"V6 



2 -'C 

1 1 A 

1 I 5c 



(2-53) 



From the parity conservation, we find only 5 coefficients are independent; namely, 
we find the correspondence 



~ /(OO) , 

[aS]^^^ + [aS]^_}l ~ 2/(11) 
[aS]f + [aS]^^l ~ 2/(22) 



(1) 



(1) 



[aS]f - [aSn ~ 



(2) 



[aS]^^^ ~ /(20) . 



(2) 



aSjf^ - [aSj'Li ~ 2/(21) 



(2) 



(2-54) 



The transformation to the Cartesian representation is carried out by using the for- 
mula 



1 

71 



(a-S) 



cr 



(2-55) 



and the formula in Eq. (|A-12p . We find the correspondence 



(0) 



(a^Sx + o-ySy + azSz 



[(75] + [aStl = [ax S]y = a^Sx - a^S, , 

\aSf'^ - {aStX = (-Ok X S]x = {-i){aySz - a.Sy) , 

^^■^^^^^ = ^ ^I*^ X 5]^ = i -^{'^xSy - aySx) , 

[aS]f^ + [aS]^^l = axSx - CFySy , [(yS]f^ - [crSf3z = {'^xSy + (JySx) , 



(2-56) 
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If we combine Eq. ()2-54p and Eq. ()2-56p and express Ca,i3 by I{Xfi), we obtain the 
following results. 

C.,. ~ -^^(00) - ^^(20) + ^(22) , 

Cy,y - -^71^(00) - ^71^(20) " ^(22) , 

C.,.~-^A00) + y|/(20) , 



C,,,~/(ll)-/(21) 



C,,, ~ -7(11) - 7(21) 



(2-57) 



Let us move to the sj^"^ = [crS'^^^]^^^ type spin correlation coefficients. The 
reduced matrix element in this case is 



(A) 



■1)1+^^/5 



1 I 

2 1 A 
1 I 5c 



(2-58) 



We again use the same notation 7(A/i) of Eq. (|2-52p with this reduced matrix element. 
Another convenient representation is decoupled representation 

= ^(la2/x|Am) Tr {/k5(2)]W/t| = ^(ia2/.|Am) 7(Am) . (2-59) 



Am 



Am 



The symmetry, = {—Y^'^^^Iatn is the same as that of Eq. ()2-45p for the 

nucleon to deuteron polarization transfer. This allows us to use the relationship in 
Eq. (|2-46p by just replacing Ka^ and Co-^ by Ca"' and 7o-^, respectively: 



C,-'^~i-^(7i2-7_i,2) , 



3 3 \ 



q.- = - {c^'- + cy'y) ~ i 2^3 7io 

_3_ 
72 

Cf^~z(-3)7o2 , 
~^37oi . 



y -'^ ~ i^n - h-i) 



(2-60) 
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We need seven independent components, Iqi, /02, ho, hi, h-i, I12, I-i,2, which 
are related to 1(20), 7(11), 7(21), 7(31), 7(22), 7(32), 7(33) through the CG coef- 
ficients. Note that 7(10) = 7(30) = 0, since 7(A0) = (-)^+°7(A0). Equation ([Ml) 
is exphcitly written as 



7l2 

1-1,2 

hi 

h,-i 
ho 

7o2 



/(33) , 



1 

1 



7(20) 



7, 



01 



7(11) + ^7(21) + a/^7(31) 
10 V2 V 5 ' ' 



?/(22) + -L/(32) 



j-/(ll)--L/(21) + 2^A,(31) 



If we use these in Eq. (j2-60p , we eventually obtain 

c^.'y ~ i 



-|7(33)-3^A7(ll) + ^^7(21) 



:/(31) 



cy/ ~ i 



Y|/(22) - ^/37(32) - ^7(20) 

7(33) + y|7(ll) - ^7(21) - y^/(31) 



3 

7^ 



--^7(33)- 2^^7(11) -3^^7(31) 



y|7(ll) + ^/67(21)+2y|7(31) 



-y|7(22)-^/37(32) + ^7(20) 



x,y 



C: 



CI' 



\/67(22) - \/37(32) 
3 



-3W— 7(11) 
10 ^ ^ 



^7(21) + 2^^7(31) 



(2-61) 



(2-62) 
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§3. Results and discussion 

3.1. Vector analyzing-power of the nucleon 

The vector analyzing-powers of the nucleon, predicted by model fss2, are shown 
in Figs. [Hand [2] for the neutron incident energies En = 3 to 65 MeV. The calculations 
in this paper were carried out using the maximum angular momentum for the NN 
system, /max = 3 or 4, and the momentum mesh-points n = 5-6-5 or 6-6-5, in the 
definition defined in I. For the energies En < 3 MeV, the partial waves up to /max = 3 
are good enough. We should note that the polarization observables are more sensitive 
to the truncation of the model space than the differential cross sections. They are 
also very sensitive to the Coulomb effect in the present low-energy region. Almost 
all polarization data are for the pd or dp scattering and the detailed comparison 
with the experiment requires the introduction of the Coulomb force. Here we show 
with dashed curves preliminary results obtained by applying the Vincent and Phatak 
methocP^ to the pd scattering^ The pd calculations were made using the cut-off 
Coulomb potential with the the cut-off radius /?c = 9 fm (for Ep < 3 MeV) or Rc = 8 
fm (for Ep > 5 MeV), together with /max = 3 and n = 6-6-5. The nd data shown 
with bars should therefore be compared with the solid curves, while the pd data 
with circles and others correspond to the dashed curves. The cited paper for these 
experimental data is specified by two letters and double figures, which are short for 
the first author's name and the publication year, respectively. In the forward angular 
region with ^cm ^ 30°, the enhancement of Ay for the pd data is almost correctly 
reproduced by the Coulomb effect, although the reduction in 9cm = 60° - 120° makes 
the agreement with the experiment worse at the energy region Ep = 3 - 14 MeV. We 
find that the difference between the solid curves and the dashed curves gradually 
diminishes for higher energies except for the forward angles. 

Our results in Fig.[T] imply that the long-standing ^y-puzzle for the large dis- 
crepancies between the theory and experiment in the energy region E^ < 25 MeV 
is not so serious as the AV18 potential,13^i3 although the maximum peak height of 
Ay is still too low. We compare in Table |T] the theoretical peak heights with the nd 
experimental data and with the pd data given in Ref. l3ip . The comparison of the 
no-Coulomb calculations implies that our results are about 80 - 90% of the observed 
nd data, although exact evaluation of the ratio is not easy due to the experimental 
errorbars. A similar amount of discrepancy is also seen in the comparison of the pd 
data, but it has an apparent minimum around Ep = 7 MeV of about 83%. As an 
average in this energy region, the shortage of about 15% is a fair estimate both for 
the nd and pd data, which is better than the rather constant discrepancy by AV18 
potential at the 25% level up to about 25 MeV.^J (See Fig. 2 of Ref.[I7j).) Looking 
back to the old predictions by realistic separable potentials. Ay (9) at some energies 
are very well reproduced.'^^^ Our QM A^A^ interaction fss2 and the separable 
potentials are both considered to have quite different off-shell properties from the 
meson-exchange potentials, that are characterized by the strong nonlocality of the 



*' The details of this approach to the pd scattering wiU be reported in a separate paper. See 
also Ref.irg)). 
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Table 1. Comparison of the maximum values of Ay and their positions Ocm with the nd and pd 
experimental data. The ratio between the theory and experiment is also shown. 



En system fss2 ratio experiment Ref. 



(MeV) 




dcm (deg) 


max 




9cm (deg) 






2 


nd 


95 


0.0341 


0.81 


95 


0.042 


EH [1.9 MeV] 




pd 


92 


0.0320 


0.973 


95 


0.0329 


EI 


2.5 


nd 


100 


0.0412 


— 










pd 


98 


0.0387 


0.921 


99 


0.0420 


ED 


3 


nd 


104 


0.0487 


0.81 


104 


0.060 


[3Sli 




pd 


103 


0.0443 


0.899 


104 


0.0493 


ED 


5 


nd 


116 


0.0726 


0.80 


115 


0.091 


[32ll 




pd 


113 


0.0684 


0.840 


113 


0.0814 


El) 


7 


nd 


118 


0.104 


— 










pd 


119 


0.0930 


0.830 


119 


0.112 


ED 


9 


nd 


122 


0.131 


0.87 


119 


0.15 


ISH) [8.5 MeV] 




pd 


123 


0.118 


0.831 


122 


0.142 


Ell 


10 


nd 


124 


0.144 


0.85 


126 


0.17 


Ml 




pd 


124 


0.131 


0.845 


124 


0.155 


EI 


12 


nd 


127 


0.167 


0.88 


128 


0.19 


[nu 




pd 


127 


0.155 


0.856 


126 


0.181 


ED 


14 


nd 


129 


0.185 


0.84 


130 


0.22 


El) [14.1 MeV] 




pd 


129 


0.177 


0.876 


128 


0.202 


El 


18 


nd 


133 


0.203 












pd 


133 


0.205 


0.928 


132 


0.221 


ED 


19 


nd 


133 


0.205 


0.89 


132 


0.23 






pd 


134 


0.208 











interaction in the configuration space. For the energies Ep = 30 - 65 MeV, there is 
no clear discrepancy of Ay between the theory and experiment. 

We should mention that the improvement of Ay (9) is not achieved by the mod- 
ification of the ^Pj phase shifts of the A^A^ interaction, as claimed in early studies 
of this problem. Although Ay (9) is very sensitive to these phase shifts, an artificial 
modification (of almost 10%) is far beyond acceptable from the modern phase shift 
analysis of the A^A^ inter action-fSJ-CB Our A^A^ model fss2 reproduces the empirical 
^Pj phase shifts within the accuracy of one degree at the energies less than 300 MeV. 
(See Fig. 1 of Ref.|T|).) From the detailed phase shift analysis of the pd scattering up 
to Ep = 10 MeV, the authors of Ref.lT5|) claim that the and £3/2- eigenphase 

shifts should be affected by the 3A^ force, in order to improve Ay (9). In fact, the 
phase shift analysis of the pd scattering is very ambiguous in this energy region, 
since so many complex parameters are involved. In Ref. ll3|) . we have compared the 
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nd phase shifts with the phase shift analysi^^ at much lower energies less than 3 
MeV and found that the desirable feature to reproduce the sufficient binding energy 
of the triton and the spin-doublet nd scattering length "^a^d is achieved by the more 
attractive feature of the "^81/2 phase shift of fss2 than the AV18 potential, originating 
from the nonlocal description of the short-range repulsion. We can conjecture that 
a similar situation is taking place even for the higher energies, although it is difficult 
to pinpoint some particular eigenphase shifts. 

3.2. Vector and tensor analyzing-powers of the deuteron 

We show the vector analyzing-powers of the deuteron in Fig.[3l and the tensor 
analyzing-powers in Figs.[l]-[6l All of these observables are measured from the dp 
elastic scattering using the polarized deuteron, and discussion including the Coulomb 
force is definitely necessary for the comparison between the theory and experiment. 
The introduction of the Coulomb force to the vector analyzing-power iTn in Fig.[3l 
shown with the dashed curves leads to the enhancement at 9cm ^ 30° and the 
reduction at ^cm = 80° - 120° for the low-energies Ep = 3 - 9 MeV, that are very 
similar to the situation of Ay. We have unpleasant rise of the hill at the angular 
region Ocm = 20° - 90° for the energies E'p = 3 - 10 MeV, which is probably related to 
the inadequacy of the nuclear-Coulomb interference term. The peak height around 
9cm = 120° is somewhat too low at Ep = 3 - 9 MeV, although it is not so serious 
as in Fig. 5 of Ref. l23p . Except for these, there is no clear discrepancies between the 
theory and experiment. As to the first problem, the enhancement by the Coulomb 
effect is also seen in other calculations by Berthold et al.,® Alt et al.,!^ Kievsky 
et al.^ and Deltuva et al.,'^''^ although to less extent in the last two cases. It is 
possible that the corresponding nd quantities are too large. 

The tensor- type analyzing-powers T20, T21 and T22 for the energies Ep = 3 
and 5 MeV in Fig. [J] show a fairly large Coulomb effect in the whole angles. In 
particular, the forward behavior of T2Q and T21 in 6cm < 60° for the no-Coulomb 
calculation (solid curves) is entirely modified by the Coulomb effect, resulting in a 
good agreement with the dp experimental data. We have some problems with T20 
and T21 in the angular region 9cm = 20°-70° for the energies Ep = 7 - 9 MeV, 
which might again related to the nuclear-Coulomb interference. In T22, too large 
negative values at the minimum points are shifted to the smaller direction, giving 
better reproduction of the experimental data. On the low-energy side with Ep < 3 
MeV, the minimum points of the dips are slightly too high. Related to this, the 
perfect fit of T22 in Ref. flU]) at E'lab < 10 MeV is fortuitous, because this calculation 
does not include the Coulomb force. The Coulomb modification diminishes smaller 
and smaller for the higher energies Ep = 7 and 9 MeV in Fig.[5l and for Ep = 22.7 
and 65 MeV in Fig.[6l The big difference from the experiment in 65 MeV T21 at 
9cm = 50° - 90° is a common feature with other calculations In order to examine 
the adequacy of our Coulomb treatment, we have compared our Coulomb effects with 
the results of calculations by other authors For example, Ref.i24j) shows 
the differential cross sections and analyzing-powers of -Ep = 3 and 9 MeV in Figs. 9 
and 10, respectively, for the coupled-channel potentials of CD Bonn, including the 
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Fig. 1. The nucleon analyzing-power Ay(Q) of the nd elastic scattering (solid curve) from En = 3 
to 19 MeV, compared with the experiment (bars). The pd results with the cut-off Coulomb 
force are also shown by dashed curves, which should be compared with empty circles. The 
experimental data are taken from Refs. [30ll for Mc94 (nrf),|3T|) for Sa94 (pd),|32]) for To91 (nd), 
ES} for Ho87 (nd),[31I for Ni98 (nd), and ESI for WelO (nd). 



22 



K. Fukukawa and Y. Fujiwara 



0.2 



0.1 



En = 22.7 MeV 

Lax=4 n=6-6-5 

* Co78 (pd) 









-0.1 



30 60 90 120 150 180 
6 cm (deg) 

Ay(e) 



0.2 
0.1 


-0.1 
-0.2 
-0.3 



E„ = 30 MeV 
\^^=4 n=5-6-5 


f \ 

\ ' 


\ i 

" 1 Do78 (nd) [29.6 MeV] \ | 
" + Jo65 (pd) [30.2 MeV] 





30 
0.4| ^ 



60 90 120 150 
6 cm (deg) 



180 



E„ = 35 MeV 




60 90 120 150 
6 cm (deg) 



180 



^(8) 
0.4 



0.2 



-0.2 



-0.4 



E„ = 40 MeV 

l^3,=4 n=5-6-5 


f K 

1 *K 


\ 

\ ' 

\ 'I 
- ^ Bu68 (pd) \ W 
\ '/ 


> 



30 60 90 120 150 180 
9 cm (deg) 

Ay(e) 

0.4r 



0.2 



-0.2 



-0.4 



En = 50 MeV 

lmax=4 n=5-6-5 




\j 1 

Ro82 (nd) \t | 

V 





30 

Ay(e) 

0.4 



60 90 120 150 
9 cm (deg) 



180 



0.2 



-0.2 



-0.4 



-0.6 



y \ 

a. 


E„ = 65 MeV A, 
\^=A n=6-6-5 /jj X - 




\ / 




\ k 


■ f Sh82(pd) 




[64.8 MeV] 





30 



180 



60 90 120 150 
9 cm (deg) 

Fig. 2. The same as Fig.[l] but for the energies En ~ 22.7 to 65 MeV. The nd data are shown by 
bars and the pd data by filled circles with bars. The experimental data are taken from Refs. l44|l 
for Co78 (pd),|13 for Do78 (nd),^ for Jo65 for Za73 {nd),\^ for Bu68 (pd),|lSll 

for Ro82 (nd), and[50j) for Sh82 (pd). 



Quark-Model Baryon-Baryon Interaction Applied to the nd Scattering (II) 23 



iTii(e) 

0.03| 



iTii(e) 



0.02 



0.01 



E„ = 3 MeV 




30 60 90 120 150 180 
d cm (deg) 

iTii(e) 



0.04 



0.02 



En = 5 MeV 




lmax=4 n=5-6-5 


3. 'Til ■ 

r \ 
4 Xl 


dp 

So87 / 
o Sh95 f 


§ % 









iTii(e 

0.08 



0.06 



60 120 
9 cm (deg) 



180 



0.04 



0.02 



- En = 9 MeV 




- \^=A n=5-6-5 


n iTii - 

/\ : 


o Sh95(dp) 


/ \° 








/ \p 

/ V 


^ '^Oaaa.^Sr — ^ 









0.05 



bn = 1 u Mev 




■ lmax=4 n=5-6-5 




1 Sp84 (dp) 




■ ^"V^-^i^ r 













I 



30 60 90 120 150 180 
9 cm (deg) 

iTii(9) 



0.2 E„ = 22.7 MeV 



5-5 




iTii 













-0.1 



30 60 90 120 150 180 
9 cm (deg) 

iTii(9) 



0.2 




30 60 90 120 150 180 30 60 90 120 150 180 

9 cm (deg) 9 cm (deg) 
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A degree of freedom. We find that the general tendency is common, although some 
quantitative differences exist because of the different treatments of the Coulomb 
force. 

3.3. Analysis of the Coulomb effect below the deuteron breakup threshold 

In our previous paper/^' we have examined the Coulomb effect on the pd differ- 
ential cross sections for the incident energies below the deuteron breakup threshold. 
In this low-energy region, the channel spin Sc is an almost good quantum number, 
and the J-averaging of the eigenphase shifts with respect to the definite Sc and the 
orbital angular momentum i between the nucleon and the deuteron is very convenient 
to reduce the number of phase shift parameters, at least in discussing the angular 
distribution of differential cross sections. However, a simple prescription adding the 
Coulomb amplitude to the nd scattering amplitude only with the Coulomb phase- 
shift factors, which is called the "Coulomb externally corrected" approximation in 
Ref. l2^ . does not work for the differential cross sections in the low-energy region. 
This prescription largely overestimates the differential cross sections. This implies 
that the modification of the nuclear phase shifts by the Coulomb force is very impor- 
tant below the deuteron breakup threshold. In Ref. ll3p . the Coulomb modification 
to the nuclear phase shifts is therefore incorporated, using the difference of the J- 
averaged eigenphase shifts for the nd and pd scatterings first from the AV18 potential 
in Ref. l22p . secondly from the fss2 calculated here. In both cases, we have obtained an 
almost complete reproduction of the differential cross sections by this prescription. 

Here, we incorporate the Coulomb force by the Vincent and Phatak method, 
directly to the nuclear scattering amplitudes. The results are shown by dashed 
curves in Figs.[7]-[9l We find that the pd or dp experimental data are generally well 
reproduced. In more detail, we find the following two problems: 

1) The Ay puzzle is more serious below the deuteron breakup threshold, especially 
for the nd data. For instance, the recent Ay measurement of the low-energy nd 
scattering in Ref.[35|) shows the discrepancy of more than 30%, although the 
difference in the pd case is not so serious. 

2) The minimum points in the deuteron tensor analyzing-power T22 are too high, 
which is a common feature with the 3 MeV result in Fig. HI 

We can compare our results with those by the variational approach in Ref. l22p . 
which incorporates the complete Coulomb force to the AV18 potential and the Ur- 
bana three-nucleon force. The AV18 potential yields particularly large discrepancy 
of more than 30% in this energy region not only for Ay but also for iTu, as seen in 
Figs. 1 and 2 in Ref.l22p. The differential cross sections and the tensor- type deuteron 
analyzing-powers are well reproduced. As to the difference in T22 in 2) above, the 
comparison with their results indicates that the three-body force might be important 
to reproduce the magnitude at the minimum points. The deuteron analyzing-powers, 
zTii and T2m, &t Ed = 3 MeV in Fig. [8] are very similar to their results in Fig. 12 of 
Ref. l22p . although the comparison with the experimental data is not easy because of 
the experimental errorbars. 
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3.4. Polarization transfer coefficients 

The nucleoli polarization transfer coefficients K'a {0) of the nd elastic scattering 
are compared with the pd data for £"„ = 10 MeV and 22.7 MeV in Fig.[TOj The dotted 
curves represent the kinematical rotation in Eq. ()2-40p when Ka^p = Sa,i3 is assumed. 
From here on, we neglect the Coulomb effect and compare the nd results directly 
with the pd data. This would be permissible since the experimental errorbars are still 
large for these polarization observables. We find satisfactory agreement between the 
theory and experiment. The vector-type nucleon to deuteron polarization transfer 
coefficients Ka (0) are compared in Fig. [11] at the same energies £"„ = 10 and 22.7 
MeV. The dotted curves again represent the kinematical rotation in Eq. ()2-44p when 
Ka,/3 = ^a,p is assumed. Various types of the nucleon to deuteron polarization 
transfer coefficients are compared with the pd data in Figs.[T2]- [TH Here, again we 
obtain satisfactory agreement although the experimental errorbars are rather large. 

3.5. Spin correlation coefficients 

The spin correlation coefficients, Cxx, Cyy, and S = (— 1/2)C^^, for the nd elastic 
scattering at En = 8.7, 10, 12 and 13 MeV are compared with the pd experimental 
data in Figs.llSI and 1161 The experimental data are actually for the dp scattering at 
the energies Ed = 17 A, 19.5, 23.8 and 26.1 MeV.El We find a reasonable agreement 
between the theory and experiment, although the experimental errorbars are rather 
large. 

§4. Summary 

One of the most important purposes of studying the three-nucleon (3A^) system 
in terms of realistic nucleon-nucleon (NN) interactions is to clarify the off-shell 
properties of the interaction, which can never be known from the physical observables 
of the two-nucleon system. It is therefore very crucial that the A^A^ interaction to 
start with can reproduce the deuteron properties and the A'^A'^ phase shifts very 
accurately. Our quark-model NN interaction fss2 satisfies this criterion by taking 
into account the naive three-quark structure of the nucleon. It is formulated in the 
framework of the resonating-group method (RGM) , in which the antisymmetrization 
of quarks generates a strong nonlocality and characteristic energy dependence of the 
interaction between two three-quark clusters. When this interaction is applied to the 
3A^ system, the short-range repulsion originating from this nonlocality behaves quite 
differently from the meson-exchange potentials in favor of the deuteron distortion 
in the spin-doublet channel. We have found that fss2 yields sufficient attraction 
comparable to the AV18 plus Urbana 3A^ force, leading to the nearly correct triton 
binding energy and the empirical value of the spin-doublet effective length for the 
nd interaction.!^ Note that this strong deuteron distortion effect is only for the 
J'^ = 1/2+ channel. On the other hand, the distortion effect is marginal in the 
spin quartet channel or the = 3/2+ channel, owing to the Pauli principle on 
the nucleon level. As the result, the eigenphase shift of the dominant ^5*3/2 state 
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Fig. 10. The nucleon polarization transler coefficients (9) of the nd elastic scattering for En = 
10 MeV and 22.7 MeV, compared with the experiment. The experimental pd data are taken 
from Ref.[53| for the 10 MeV data and from Ref.|61]) for the 22.7 MeV data. The dashed bars 
in the panel of 22.7 MeV Ky^' are from Refl52]). 
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Fig. 11. The vector- type nucleon to deuteron polarization transfer coefficients (9) of the nd 
elastic scattering for En — 10 MeV and 22.7 MeV, compared with the pd experimental data. 
The experimental data are taken from Ref.l53[l for the 10 MeV data and from Ref.l63p for the 
22.7 MeV data. 
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Fig. 12. The tensor-type nucleon to deuteron polarization transfer coefficients ~' {6) of the nd 
elastic scattering for En — 10 MeV and 22.7 MeV, compared with the pd experimental data. 
The experimental data are taken from Ref.l53p for the 10 MeV data and from Refs.l61 [) .l63 [) for 
the 22.7 MeV data. 



Quark-Model Baryon-Baryon Interaction Applied to the nd Scattering (II) 35 




36 



K. Fukukawa and Y. Fujiwara 



predicted by fss2 is very similar to the AV18 potential, in which the effect of 3A^ force 
is very small. In the low-energy region, the effect of the 3A^ force hardly appears since 
the differential cross sections are dominated by the partial waves of the spin-quartet 
channel. It is important to note that fss2 can also reproduce the differential cross 
sections at higher energies, including the sufficient magnitude of the cross sections 
at the diffraction minima.^-' 

In this paper, we have extended the previous studies,'^''^ applying fss2 to the 
nd elastic scattering, to various types of polarization observables. The long-standing 
Ay-puzzl^^^^ for the nucleon vector analyzing-powers in the low-energy region 
En ^ 25 MeV is largely improved in comparison with the predictions by the AV18 
potential. Although the Coulomb effect obscures the definite conclusion, the short- 
age of the maximum peak is about 15%, which is somewhat similar to the predictions 
by old realistic separable potentials .^D^® The detailed comparison with the phase 
shift analysis is not possible since the empirical phase shifts are not uniquely deter- 
mined. Nevertheless, the conclusion drawn in Ref. ll6p based on the solutions derived 
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Fig. 14. The same as Fig. ll2l but for the other types of tensor-type nucleon to deuteron polarization 
transfer coefficients. 
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from the AV18 potential should be reconsidered, since ^Pj phase shifts at the energy 
Ep ~ 10 MeV are quite different between fss2 and the AV18 potential. The difference 
is of the order of several degrees. It is possible that this difference is caused by the 
nonlocality of the quark-model A^A^ interaction. Since all the experimental data for 
the deuteron analyzing-powers are for the pd or dp scattering, we have introduced the 
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cut-off Coufomb force by the Vincent and Phatak method. This method works very 
weh to reproduce the behavior of the differential cross sections and analyzing-powers 
at the forward angles. In particular, the low-energy observables below the deuteron 
breakup threshold are well reproduced except for the ^y-puzzle and a slight under- 
estimation of the dips in T22{0) at the minimum points around ^cm = 100° - 120°. A 
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Fig. 16. The same as Fig.[TSl but for E„ = 12 MeV and 13 MeV. 
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large discrepancy of the vector- type analyzing-power of the deuteron, iTii{6), found 
in Ref. [23|) is not observed in our calculations, although the peak height around 
Scia = 120° is somewhat too low at E'p = 3 - 9 MeV. Instead, we have unpleasant rise 
of the hill at the angular region ^cm = 20° - 90° for the energies Ep = 5 - 10 MeV, 
which is probably related to the inadequacy of the nuclear-Coulomb interference 
term in our Coulomb treatment. A similar problem is also seen in T20 and T21 in the 
angular rigion ^cm = 20° - 70° for the energies Ep = 7 - 9 MeV. Except for these, 
there is no clear discrepancies between the theory and experiment. We have care- 
fully examined the Coulomb effect on the observables between our results and other 
calculations using the modern meson exchange potentials like AV18 and CD Bonn 
potentials .l23~ESli fpj^g direction of the Coulomb modification is always the same, al- 
though some quantitative diference appears at the forward angles. Our calculation 
somehow overestimates the Coulomb effect and nuclear-Coulomb interference term 
is not precisely reproduced. On the other hand, the behavior of the observables in 
the backward angles 6cm ^ 90° is hardly influenced especially at higher energies 
Ep > 10 MeV. We have also examined various types of polarization transfer coeffi- 
cients and the spin correlation coefficients by neglecting the Coulomb force. They 
are reasonably reproduced, although the experimental errorbars are still very large. 

In conclusion, we find no apparent disagreement between the theory and exper- 
iment as far as nd elastic scattering with the energies En < 65 MeV is concerned. 
It should be stressed that this conclusion is valid only when we deal with the en- 
ergy dependence of the quark-model RGM kernel properly.!^ As discussed in 
our previous paper,® this energy dependence is eliminated by the standard off-shell 
transformation utilizing the square root of the normalization kernel. This proce- 
dure yields an extra nonlocal kernel which is not extremely small and affects various 
3A^ observables in different ways. In the next paper, we will discuss the deuteron 
breakup pro cesses .fSZli 
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Appendix A 

Spin operators of the deuteron 



In this appendix, we summarize various notations for the spin operators of the 
deuteron. The deuteron spin operators are most easily defined from the Wigner- 
Eckart theorem 



(A-1) 
VT^. Here, 



with the standard value of the reduce matrix element (l||S||l)un( 
the subscript "unc" stands for the unconventional reduced matrix element with 
(S'||)S||S')unc = ^/S{S'^\-l). This gives the rank-one spin operator of the deuteron as 

(A-2) 

The transformation to the Cartesian representation is carried out by the standard 
spherical vector representation: 
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(A-3) 
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(A-4) 



Note that 5"^ = 5 and 5^ 

The rank- two spin operator of the deuteron is defined by Sm^ = [SS]l 



2, which is consistent with the reduced matrix element. 

:(2) _ rcc](^) and the 

reduced matrix element is calculated from 



(A-5) 



(W|5(2)|lm) = {lm2fi\lm'){l\\S^^^\lU,, . 

We note that the spin operator S does not change the spin value S = 1, when it is 
operated on the deuteron spin wave function xi = Xi(12): 

S\Xl)(ll)l^^ = -[SxiW = ^XiM • (A-6) 
(Note the phase change by the angular-momentum coupling order.) We find 

"112 

s^'^xi)ii2)i, = [s^'hiW = [[sshxiW= E oil 

A=0,l,2 J ^ J 



[S[Sx 
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1 

1 1 



2 1 
1 1 



(_V2)(-^/2)xi^ 



which yields 

This value and Eq. ()A-5p gives 

(W|5(2)|lm) = 



1— m 



{Im'l — m|2/i) , 



(A-7) 

(A-8) 
(A-9) 



owing to the symmetry property of the CG coefficients. The explicit values of the 
CG coefficients 



(1111|22) = 1 , (1110|21) 
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and other non-zero coefficients obtained by the symmetries, yield 
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The transformation to the Cartesian representation is carried out by using a general 
formula which is obtained from the CG coefficients Eq. (jA-lOp and the vector-tensor 
transformation Eq. (|A-3p . Namely, for two vectors a and b we have 

[ab]22 + [ab]2-2 = a^bx - ayby [ab]22 - [a&]2-2 = iiaxby 

[ab]2i - [ab]2-i -- 



(^ybx) 1 



[ab]2i + \ab]2-i = - 
[a6]2o = (20262 



-i{aybz + a^by) 
Oixbx 



ayby) 



[oxbz + Ozbx) , 

(A-12) 



If we set a = b = S, we find 
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Here we should note that the deuteron spin operator S satisfies the commutation 
relations 

[Sa, 5/3] = i ^ £al3-/S.y . (A-14) 
7 

We use this relationship and define the spin operators of the deuteron in the Cartesian 
representation Vai3, following the definition by Ohlsen (see Eq. (2.17) of Ref.[29]).) 
Namely, we define Va = Sa and 

3 3 

"Pxy — 2 {SxSy + SySx) , "Pxz — 2 {SxSz + SzSx) , 

ryZ = ^{SyS,+S,Sy) . (A- 15) 

For the diagonal part Vaa, we define Vaa = 35q,^ — 2 symmetrically, which leads to 
'Pxx + T^yy "l~ ~ Oi owing to 5^ = 2. We also use the normalized tensor operators 
of the deuteron Ti^ and defined by 



which satisfy 



Tj = {-YT^-u , ^Tr (t^uTI^,) = 8^,,>5,y . (A-17) 



The Glockle et al.'s review article^ also uses the notation S^^ (i = - 8) for T^^, 
which satisfy 

S^^^ = S^'^ , ^Tr [s'^^ = 6ij . (A-18) 

All of these different notations are related by 

5(0) = 1 , 



s^'^ = Jlv. = Jls. = JU 
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0-1 
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= a/I?'.. = {s^s, + s^s,) = I I -1 



0-10 



sm = -L(P„-P„) = -L3(s/-V) = y|(o 



0-10 

1 
^10 




1 
0-2 
1 



or 



(A-19) 



The relationship between the tensor notation and S^*) is 

5f ) = -L (^3^7) + -5(4)^ ^ ^(2) ^ _^ (^^(7) _ .^(4)^ ^ 

5f = , (A.20) 



Too = S('^ = 1 , 
Tio = S^^^ , 

T20 = . (A-21) 
If we express these by Va and Vap, we obtain 
1 x/s 1 

^2 (Til + = —Vy , r2o = -^P.. , 

^ (r2i - r2,_i) = -^^xz , ^ (T22 + r2,_2) = ^ {v,, - Vyy) , (a-22) 

where the other components of the deuteron analyzing-powers are zero as in § 2.2. 
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